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Abstract 

In spite of extensive observations and numerous theoretical studies in the 
past decades several key questions related with Gamma-Ray Bursts (GRB) 
emission mechanisms are still to be answered. Precise detection of the GRB 
polarization carried out by dedicated instruments can provide new data and 
be an ultimate tool to unveil their real nature. A novel space-borne Compton 
polarimeter POLAR onboard the Chinese space station TG2 is designed to 
measure linear polarization of gamma-rays arriving from GRB prompt emis¬ 
sions. POLAR uses plastics scintillator bars (PS) as gamma-ray detectors 
and multi-anode photomultipliers (MAPMTs) for readout of the scintillation 
light. Inherent properties of such detection systems are crosstalk and non¬ 
uniformity. The crosstalk smears recorded energy over multiple channels 
making both non-uniformity corrections and energy calibration more diffi- 
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cult. Rigorous extraction of polarization observables requires to take such 
effects properly into account. We studied influence of the crosstalk on energy 
depositions during laboratory measurements with X-ray beams. A relation 
between genuine and recorded energy was deduced using an introduced model 
of data analysis. It postulates that both the crosstalk and non-uniformities 
can be described with a single matrix obtained in calibrations with mono- 
energetic X- and gamma-rays. Necessary corrections are introduced using 
matrix based equations allowing for proper evaluation of the measured GRB 
spectra. Validity of the method was established during dedicated experimen¬ 
tal tests. The same approach can be also applied in space utilising POLAR 
internal calibration sources. The introduced model is general and with some 
adjustments well suitable for data analysis from other MAPMT-based in¬ 
struments. 

Keywords: Muti-anode photomultiplier; Non-uniformity; Crosstalk; 
Gamma-ray burst; Polarization; 


1. Introduction 


Gamma-ray bursts (GRBs) are observed as short flashes of gamma-rays 
appearing randomly in the sky. In a few seconds they release energy be¬ 
tween J and 10^® J making them the most energetic explosions in the 
Universe. They are frequently associate with either a collapse of the mas¬ 
sive stars or a violent merge of compact binaries [H. In spite of numerous 
observations and theoretical efforts in the past decades many key questions 
such as GRB emission mechanisms or origin and structure of their magnetic 
held are not answered yet (for recent reviews, see Refs. e.g. 0 , 31 ) . Direct 
measurements of the GRB polarization in the prompt emission phase should 
be able to shed light on the whole system and constrain the energy emission 
mechanisms 0, IJ 

To achieve this goal, several dedicated instruments are currently under 
development 0,3 3 3 113 • They are specially designed and optimized for 
polarization measurement with support from precise on-ground calibrations 
as well as rigorous performance modeling and verihcation. POLAR is one of 
such new hard X-ray polarimeters aimed to study the GRB prompt emissions. 
It utilizes Compton scattering to measure linear polarization of gamma-rays 
in the energy range from 50 keV to 500 keV. Plastic scintillator (PS) bars 
are chosen as gamma-ray detecting medium and 8x8 channels multi-anode 
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photomultipliers (MAPMTs) serve for readout of the scintillation photons. 

It is an inherent property of MAPMT-based detectors that a signal from 
a certain channel can, through crosstalk effects induce signals in the neigh¬ 
boring channels. In addition, the response non-uniformity of the MAPMT 
modifies initially deposited energy altering the signal amplitudes. Together 
with the crosstalk effect that spreads the initial energy deposition, it makes 
the energy calibration and spectral unfolding more difficult. Thus, the pre¬ 
cise knowledge of both effects is necessary to properly extract polarization 
observables in the detected GRBs. 

Based on laboratory calibration data we constructed a model describ¬ 
ing both the crosstalk and response non-uniformities for all 64 channels of 
each POLAR module. The relation between genuine and recorded energy 
deposition can be described by a single matrix. In the following chapters we 
describe methods applied to determine matrix elements and present results 
from laboratory tests used for their verification. 

2. POLAR instrument 
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Fig. 1: Exploded view a POLAR detector module (left) and the full POLAR instrument 
(right). Each module has 64 PS bars (5.9 x 5.9 x 176 mm^ each) coupled to a 64 channel 
MAPMT (Hamamatsu H8500) and a front-end electronics. The full instrument consists 
of 25 identical modules. 

The main goal of POLAR to measure linear photon polarization is real¬ 
ized using Compton scattering. Polarized gamma-rays undergoing Compton 
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process tend to scatter perpendicularly to their incident polarization vector 
according to the Klein-Nishina equation: 


da 

dil 




^ — 2 sin 9'^ cos 
E 


( 1 ) 


where is the classical radius of the electron, E and E' are the energy of 
the incident photon and the scattered photon, respectively, 6 is the scatter¬ 
ing angle between initial and hnal photon direction, and r] is the azimuthal 
scattering angle between the initial polarization vector and the direction of 
the scattered photon. 

POLAR has both, a large effective detection area (~ 80 cm^) and a wide 
held of view (~ 1/3 of full sky). They are needed for efhcient and precise mea¬ 
surements of the azimuthal distribution of gamma-rays undergoing Compton 
scattering in its scintillator bars. After being hit by a gamma-ray the instru¬ 
ment records energy depositions in its 1600 channels. POLAR uses an array 
of 40 (row) X 40 (column) plastic scintillator (PS) bars as gamma-ray de¬ 
tection target. The azimuth angle of the scattered gamma-ray is determined 
from positions of two bars with the highest energy depositions. Polarization 
degree as well as polarization angle of the detected GRB can be retrieved in 
the off-line reconstruction of all recorded gamma-ray events. 

The scintillating material EJ-248 was chosen because of its fast response 
and high value of the softening temperature (90°C). Each PS bar has dimen¬ 
sions of 5.9 X 5.9 X 176 mm^. In order to reduce optical crosstalk all bars 
have narrower bottom-end cross-sections resembling a pyramid-like shape. 
The procedure used for bar selection involved several strict criteria. Firstly, 
all of them were carefully inspected to reject macroscopic defects. Secondly, 
the dimension of each bar was precisely measured and only bars with dimen¬ 
sion deviations smaller than 0.1 mm were accepted. Afterwards, the light 
output difference between the top and the bottom of each bar was measured 
with a dedicated setup. It consisted of a photomultiplier and an Aluminum- 
made bar holder lined with the Enhanced Specular Reflector (ESR) Elms. 
The test used a collimated Am-241 source placed at the bottom and the 
top of the bar. Only bars with the light output difference smaller than 10% 
were accepted. Finally, to increase the light collection, all selected bars were 
wrapped in the 65 pm thick ESR him. 

The 1600 selected bars were assigned to 25 identical modules. Each mod¬ 
ule consists of 8x8 PS bars, a soft optical coupling pad made of transparent 
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Fig. 2: Numbering convention of POLAR module channels adopted for data analysis from 
Hamamatsu H8500 MAPMT [HI (seen from the top of the glass window). The dimensions 
of the MAPMT channels are also shown. 

silicon (Dow Corning DC93-500), a 64 channel flat panel MAPMT (Hama¬ 
matsu H8500) and a front-end electronics (FEE). This structure was en¬ 
closed in a 1 mm thick carbon fiber sockets (see the left pannel of Fig. [T]). 
The top and the bottom of each PS were fixed and aligned with two plastic 
frames. They enhance resistance to vibrations and further reduce the optical 
crosstalk. The dimension of the optical coupling pad is 50 x 50 x 0.5 mm^. It 
can partially absorb vibrations and also protects the MAPMT glass window. 
The channel numbering convention and MAPMT dimensioning are shown in 

Fig. [3 

The front-end electronics consists of three stacked Printed Circuit Boards 
(PCBs). The boards contain low voltage power supply circuits, a dedicated 
voltage divider for the MAPMT, a multi-channel ASIC (IDEAS VA64), an 
ADC, a FPGA and a temperature sensor. A special internal pulser circuit is 
also included to test the throughput and gain of each input channel. More¬ 
over, at the inner edge of each corner-located module a low activity, point-like 
22-Na radioactive source is installed. Two annihilation photons (511 keV) 
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from the source are used for the in-flight calibration done by applying the 
same data analysis techniques as outlined below. If the incoming signal has 
an amplitude above the threshold the ASIC sets a trigger flag that is sensed 
by the Central Task Processing Unit (CT). If CT accepts the event, the 
front-end electronics digitises previously hold amplitudes from the 64 analog 
signals, packs them into the Science Data Packet and transmits to CT for 
further processing. In addition, the arriving time and the hit pattern of the 
trigger flags are also recorded by CT and form the Trigger Data Packet. Ded¬ 
icated part of CT also manages the low and high voltage power supply and 
the settings for all FEEs. It also handles communication with the spacelab. 
All 25 modules, the CT, the power supplies and interfacing electronics are 
enclosed in a housing box as shown in the right panel of Fig. [T] The whole 
instrument is mounted on the outside pannel of of the Chinese TG-2 spacelab 
facing to the outer space. 

Monte Carlo simulations showed that for each Compton scattered photon 
the line connecting two PS bars with the maximum energy depositions is 
correlated with the outgoing azimuthal photon direction. We define ^ as the 
azimuthal angle between the line connecting these two bars and the x-axis 
of the POLAR detector. According to previous results (see e.g. Refs. 


and (l2|), the distribution of ^ from all detected GRB gamma-rays is called 


a modulation curve and follows characteristic pattern: 

/(p) = a: - {1 +/icos[2(^ -^o)]}, 


( 2 ) 


where K is the normalization factor, /i is the modulation factor and is 
the polarization angle. The polarization degree p is equal to p = p/pioo, 
where pioo the instrument modulation factor for 100% polarized gamma-rays. 
Realistic Monte Carlo simulations and calibrations with 100% polarized x- 
rays are used to determine the value of pioo- The modulation factor p and the 
polarization angle .^o are obtained from the ht to the experimental data. For 
detailed description of the reconstruction of both polarization observables 
see Refs. 


3. Crosstalk and response non-uniformity modeling 

Optical photos in PS bars are produced by ionizing particles as result 
of their energy depositions. The average number of optical photons Abar 
collected at the bottom of the bar is given by 

A^bar S • C • F/(jep, (3) 
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where i?dep is the real energy deposition (i.e. visible energy deposition), s is 
the scintillation efficiency, i.e. the average number of optical photons pro¬ 
duced per unit of deposited energy, and c is the photon collection efficiency, 
i.e. the fraction of the scintillation photons reaching the bottom of the bar. 
The cumulative coefficient c parametrizes effects of imperfections in micro¬ 
cracks, reflectivity and attenuation of the optical light in the bars. Knowing 
that all POLAR PS bars came from the same production batch and the in¬ 
flight temperature differences between them are minor, it is reasonable to 
assume that they all have the same scintillation efficiency. Hence, for each 
POLAR detector module, the number of photons reaching the bottom of its 
64 PS bars A'bar = (Abaro, Abari, ''' , cau be described as 


Abar — B.Edep, (4) 

where the vector E^ep = {Edep,o, Adep,i, • • •, Adep, 63 )'’" represents 64 real energy 
depositions in the bars, and B is a diagonal matrix expressed as 


B = sDiag(ci,C2,---,C64). (5) 

Due to crosstalk effects any signal appearing in a certain channel can 
simultaneously induce signals in its neighboring channels. According to the 
H8500 MAPMT datasheet 0. a typical crosstalk factor between two neigh¬ 
boring channels measured with an optical fiber of 1 mm diameter is of the 
order of 1%. Measurements in Ref. [l^ also supports this value. As this 
factor includes both optical and electrical contribution e.g. from stray or 
secondary electrons in the dynodes, one can conclude that the electrical part 
is smaller than 1%. On the other hand, measurements with fully assem¬ 
bled POLAR modules showed that the total crosstalk between neighboring 
channels can be as high as 10% - 20%. This contains further contributions 
from stray photons in the optical pad and spreading of photons between PS 
bars. Additional tests showed that the electric crosstalk in the module read¬ 
out system is negligible [l^. As the overall electrical crosstalk is very small 
compared to the optical one, we only consider the later. In order to describe 
its strength, we introduce a 64 x 64 matrix X = (x^), in which the matrix 
element Xij represents the fraction of photons in the j-th channel coming 
from the primary energy deposition in the i-th channel. Obviously, we have 

63 

0 < Xij < 1 and — 1- According to the principle of reversibility, we 

3=0 


7 


also have Xtj = Xji. The number of optical photons iVpm reaching all 64 
photocathodes of the MAPMT can therefore be given by 

^pm = XiVbar- (6) 

The MAPMT transforms optical photons absorbed by its photocathode 
into electric signals on its anodes. These signals are subsequently read out by 
the FEE ASIC. For gamma-ray energies detected by POLAR one can assume 
that both the MAPMT and the FEE are in their linear range of responses. 
Hence, the relationship between the recorded signal Emeas, i-e. the recorded 
energy depositions and the number of the optical photons reaching MAPMT 
A^pm can be given by 

Eraeas = GLVpjn, (7) 

where G = Diag (^fi, 5 ^ 2 , • • •, ( 764 ) represents the averaged signal induced by 

the optical photon in each of 64 channels. It is worth stressing that G 
is an individual feature of the MAPMT and FEE subsystem. In addition 
to the collection efficiency, G also describes the module uniformity level. 
Since both temperature and high voltage can change the system gain, careful 
experimental calibrations are required to obtain correct value of G. From 
Eqs. (jl]), ([6]) and ([7]), we have 

Eraeas = (GXB)Edep = R-.Edep, (8) 

where R = (rjj) = GXB is called the response matrix given by 

Tij = SQiXijCj. (9) 


Since Xij represents the optical crosstalk between two channels, both Xij and 
Tij decrease with the distance between them. 

The real energy E^ep deposited by gamma-rays in all 64 PS bars can be 
reconstructed by a linear transformation of the recorded energy depositions: 





( 10 ) 


This transformation allows for a simultaneous correction of both crosstalk 
effects and non-uniformities. By knowing the values of both E^ep and corre¬ 
sponding Ejneas oue cau also determine the response matrix R . The method 
used for its construction is described in the next chapter. It should be men¬ 
tioned that this technique can be applied either for laboratory tests on ground 
or during real observations in space. 


4. Determination of response matrix 


The energy deposited in PS bar in the process of Compton scattering is 
a function of the photon scattering angle. It leads to a spectrum of energies 
with a sharp cutoff at the end of the spectrum called the Compton edge. The 
edge is related to the maximum energy transfered to the electron and full 
back-scattering of the gamma-ray. It is commonly used for energy calibration 
in detector systems based on plastics scintillators such as POLAR. 

Let us consider a scenario in which a gamma-ray is fully back-scattered 
in the Ath POLAR PS bar and escapes from the module. The energy is 
deposited in only one i-th channel and is equal to the value of the Compton 
edge position E^e = E^\i — 1/{1 + 2i?.y/511)], where E^ is the energy of the 
incident gamma-ray in units of keV. According to Eqs. ([8]) and (jO]), the 
corresponding recorded energy deposition iff®, is given by iff® = ru ■ Ece- As 


described in Ref. [13|], iff® can be obtained by htting the recorded energy 


spectrum with a smeared step-like function. Hence, the diagonal element rj, 
of the response matrix is given by: 


r 


ii — 



( 11 ) 


The physical meaning of is the energy conversion factor of the i-th channel 
(in units of ADC channel/keV). Applying the same procedure to all other 
channels provides the energy conversion matrix M described as follows: 


M = Diag (ri,i, r 2 , 2 , • • •, ^64, 64 ) • (12) 

It should be noted that especially at lower gamma-ray energies the full energy 
absorption peaks from the photoelectric effect are frequently used for energy 
calibration of plastic scintillators. In case of POLAR, calibrations with the 
Compton edge are easier than extracting of photo-peaks due to much lower 
probability of the photoelectric process and poor PS energy resolution. In 
order to characterize POLAR response near its low energy threshold, several 
test measurements with the use of photopeaks were also conducted. 

As described above, the diagonal elements of the response matrix can be 
determined by energy calibrations. In order to determine the non-diagonal 
elements of the response matrix, let us again consider the scenario that the 
gamma-ray deposits energy only in the i-th channel. In such case, we have 
the corresponding energy deposition vector E^ep = (0, • • •, Ef^^, • • •, O)"*". The 
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recorded energy depositions of the i-th channel and of the j-th channel 

Ejieas^ according to Eqs. ([H]) and (E]), are given by 

Er^ = P, (13) 

and 

^meas _ rjiEf^^. (14) 

Since the energy deposition only occurs only in the i-th channel, the signal 
recorded in the j-th channel (i 7 ^ j) conies from the crosstalk. According to 
Eqs. dl]) and (na, the ratio fij between the recorded energy depositions of 
the i-th and the j-th channel , is called the crosstalk factor and is given by 


fij — 


^meas 

^meas 


r - • 






(15) 


The crosstalk factor fij represents partial transmission of the energy depo¬ 
sition from the i-th to the j-th channel. Note that fij differs from the pure 
optical crosstalk factor Xij as it also includes module non-uniformities. Ob¬ 
viously, we have fij = 1 for i = j, and fij 7 ^ fji for i 7 ^ j in most cases. 
One way to obtain the crosstalk factor fij is to determine the distribution 
of the recorded energy depositions of the j-th channel as the function of the 
recorded energy deposition of the i-th channel in Ref. 13|. For this pur¬ 
pose one can use pencil-like gamma-ray beams e.g. from the synchrotron 
light sources hitting exactly one PS bar. Either Compton effect or photo¬ 
absorption process can be used for crosstalk analysis. POLAR instrument 
was characterized both ways using several gamma-ray energies during test 
runs at the European Synchrotron Radiation Facility (ESRF) in Grenoble. 
It can be shown that the fit of the distribution with a line is usually suffi¬ 
ciently accurate for proper determination of the crosstalk factor. Applying 
this procedure to all pairs of channels allows to construct the full 64 x 64 
matrix F = [fij) called the crosstalk matrix. From Eqs. (jS]), (jD]), flT^ and 
flT5|) we find that the response matrix R can be calculated from these two 
matrices using the equation: 


R = F^M. 


(16) 


5. Model verification 

We have conducted a series of extensive calibration runs of POLAR at 
ESRF with synchrotron hard X-ray beams and at PSI using laboratory 
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Fig. 3: Typical energy spectra from a single channel illuminated by X-ray beams with 
energies of 356 keV and 511 keV. The Compton edges provide energy conversion factors 
of 32.1 ADC channel/keV and 33.4 ADC channel/keV respectively. The difference is 
attributed to small non-linearities of the readout electronics. 


gamma-ray radioactive sources and lower energy X-rays from the X-ray gen¬ 
erator. Some of the low energy calibration runs and the studies of the photo¬ 
absorption peaks are described in Ref. 1^. They were primary applied to 
determine responses near the low energy threshold and correct the ionization 


quenching effects that are significant at low energies 17|]. For this analysis 


we use the data from calibration runs taken with the POLAR qualihcation 
model (QM) at the ESRF ID 15A beamline with the X-ray energies of 356 keV 
and 511 keV. The synchrotron beam had a square shape of 0.5x0.5 mm^. In 
order to avoid pile-ups the beam intensity was decreased to ~ 2000 photons 
s“^ using an Aluminum wedge. One module of the QM was mounted on the 
X-Y table driven by two step motors. They allowed to change the module 
position remotely and bring the beam to the center of each PS bar. During 
the tests, each bar was illuminated for 10 seconds and then the next one was 
positioned at the beam. 

As PS has better energy resolution at higher energies, both the 356 keV 
and the 511 keV X-rays were selected to determine the energy conversion 
factors. The procedure was used to verify and correct for possible response 
non-linearities e.g. in the electronic readout system. Fig. [3] shows typical 
energy spectra recorded by a single channel when its bar was illuminated 
with photon beams at energies of 511 keV and 356 keV. The environmental 
background during the tests was negligible. The Compton edge position in 
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Fig. 4: Energy conversion factors (i.e. diagonal elements of the energy conversion matrix 
M) of one POLAR module obtained from the Compton edge fits of two energies. The 
hgure also indicates the response non-uniformity within a module. 


each spectrum was found by fitting the right edge with the step-like function: 


/(x) = oo + ai • Erfc[(x - 02 ) • 03 ], (17) 

where Erfc(a;) is the error function given by 

2 

Erfc(x) = —= 

VTT 



The parameter 02 is approximately eq ual to the position of the Compton 
edge in the units of ADC channel [13|. The energy conversion factor was 
calculated as a^lE^e-, where E^e is the Compton edge position in units of 
keV. Fig. 0] shows the energy conversion factors for all 64 channels from the 
Compton edge hts. The hgure presents the level of module non-uniformities 
and indicates some minor non-linearities in the energy response. The non¬ 
linearity effect was due to increased walking of the triggers caused by high 
values of the thresholds. It affected the instant of time of holding/sampling 
in the readout electronics. The mean value of two energy conversion factors 
was used to construct the energy conversion matrix M. 

To determine the crosstalk we applied a method similar to one described 


in Ref. ISj. In Fig. |5]we show an example of the mean energy depositions 
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(ADC channel) 

Fig. 5: Crosstalk data as mean energy depositions in channel 11 as a function of the 
energy depositions in channel 10 that was illuminated by X-ray beams with energies of 
356 keV and 511 keV. Two crosstalk factors from the linear fits at two energies are equal 
to 12%. 


recorded by channel 11 as a fnnction of energy deposition recorded by chan¬ 
nel 10. Only channel 10 was illnminated using either the 356 keV or the 511 
keV X-ray beam. Recorded energy depositions in all crosstalk channels are 
predominantly much smaller than depositions in the primary channel. Thus, 
only events for which the directly illuminated channel had larger recorded 
energy depositions were selected for data shown in Fig. |5l According to the 
Monte Carlo simulations both the 511 keV and the 356 keV X-rays have less 
than 1% chance to deposit energies larger than 5 keV in one of its closest 
neighboring bars through direct scattering processes. Thus, the contribution 
of such scattered X-rays to the crosstalk distribution showed in Fig. Elis negli¬ 
gible. The data points can be htted with a line and its slope value represents 
the crosstalk factor. The slope values for the crosstalk from channel 10 to 
channel 11 are for both energies equal to 12%. It shows that it is sufficient to 
use just one X-ray energy to determine the crosstalk factor values. Applying 
the same routine to correlation plots between any two channels of the same 
module provided the full crosstalk matrix F as shown in Fig. |61 The groups 
of channels observed in Fig. [6] are due to the fact that the channel numbers of 
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neighboring bars are not always continuous (see Fig. |2]). As shown in Fig. [3 
the crosstalk factors strongly decrease with increase of the distance between 
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Fig. 6: An example of the crosstalk matrix F of the selected module. Matrix elements 
represent crosstalk factors between two channels. Observed groups are due to the fact that 
the channel numbers between neighboring bars are not always continuous (see Fig. Ill). The 
MAPMT channel numbering convention is shown in Fig. [2j 



Distance between bars (Nb. of bars) 


Fig. 7: Mean crosstalk factor as a function of the distance between two bars calculated 
for all pairs within one module. The mean value for the closest channels is on the level of 
15%. 
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two channels. The mean crosstalk factor value between two closest neighbor¬ 
ing channels is at the level of 15%. It should also be noted that the crosstalk 
factors can be determined from measurements with either gamma-rays from 
radioactive sources or even an environmental background photons by apply¬ 
ing exactly the same method as above. For this purpose four weak ^^Na 
radioactive sources were placed inside of the POLAR Flight Model (FM). 
As they emit two collinear annihilation photons the event selection will be 
unambiguous allowing for proper calibration of the whole instrument. It will 
be possible to calibrate energies as well as determine crosstalk factors and 
non-uniformities during the entire flight onboard TG2. 
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Fig. 8: Module response matrix R calculated with Eq. (ITHll . 

Fig. [8]shows the response matrix R calculated with Eq. ffT6|) . Real energy 
depositions for two test runs were reconstructed with the help of R according 
to Eq. (HUD. Both the crosstalk and the non-uniformities were to very large 
degree corrected. As an example, the left panel of Fig. |9] shows the values of 
the residual crosstalk from channel 10 to channel 11 in the corrected data. 
We again htted residual crosstalk plots with the first order polynomial to 
determine the remaing values. The distribution of the residual crosstalk 
factors between two closest neighboring channels is shown in the right panel 
of Fig. [9l The mean value is equal to 0.01% with a standard deviation of 
0.5%. 

In order to determine the correction of the module non-uniformities, we 
selected hits produced by the beams in the PS bars. The same step-like 
function described in Eq. f[T7|) was used to ht the Compton edge positions in 
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Fig. 9: Residual crosstalk from channel 10 to channel 11 after applying corrections 
(left) and distribution of the residual crosstalk factors (i.e. the slopes of linear fits of 
the residual crosstalk plots) between two closest channels (right). The mean value of the 
residual crosstalk factors is equal to 0.01% with a standard deviation of 0.5%. The data 
from both runs at energies of 511 keV and 356 keV was used for the residual crosstalk 
studies. 


the energy spectra. Compton edge positions of all 64 channels for runs at X- 
ray energies of 356 keV and 511 keV are shown in Fig. [101 The mean values 
of the corresponding Compton edge positions in the reconstructed energy 
spectra are 208.0 keV and 341.0 keV, respectively. The theoretical values of 
Compton edge energies are equal to 207.3 keV and 340.7 keV respectively. 
From the initial maximum non-uniformity ratio reaching the value of more 
than 200% the hnal level found after correcting the data is equal to 1.2% for 
the 511 keV run and 2.0% for the 356 keV run. 
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Fig. 10; Compton edge positions of the 64 channels in the reconstructed energy spectra for 
runs at energies of 356 keV (left) and 511 keV (right). The crosstalk and non-uniformity 
corrections were applied to the selected events. Theoretical values of Compton edge ener¬ 
gies are presented as dashed lines. 


6. Conclusion 

Polarization measurements in the new POLAR instrument are based on 
the Compton scattering and utilize distribution of the azimuthal angle for 
hard X-rays scattered between its 1600 plastic scintillator bars. As the initial 
values of the POLAR response non-uniformity may reach more than a factor 
of two and the crosstalk factors could be higher than 20% one must take 
these effects properly into account. Therefore, for precise determination of 
the polarization observables in GRBs, a model describing the influence of the 
crosstalk and non-uniformities was developed. Based on it, both factors can 
be corrected for in a rigorous way. The procedure uses a linear transforma¬ 
tion of the recorded energy depositions with a specially constructed response 
matrix. A set of dedicated calibration runs was conducted for experimen¬ 
tal determination and verification of the model parameters. The procedure 
involved two steps in which the energy calibration and crosstalk measure¬ 
ments were performed using mono-energetic, pencil-like X-ray beams. The 
real energy depositions were obtained for each POLAR channel by applying 
the transformation formula as described above. To validate the method, the 
data from two test runs performed at ESRF using photons with energies of 
356 keV and 511 keV were analysed. Initial energy spectra were unfolded 
using previously constructed response matrix with implemented correction 
factors. The reconstructed data showed that the residual crosstalk factor be¬ 
tween two neighboring channels was canceled almost to zero with a standard 
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deviation of 0.5%. The remaining non-uniformity values were at the level 
of less than 2%. The presented method provides a useful tool for more reli¬ 
able computation of the real energies deposited in PS bars by the incoming 
photons. It allows for precise determination of the polarization observables 
in the detected GRB events during POLAR operation in space. With some 
adjustments this method can also be used for other MAPMT-based instru¬ 
ments. 
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